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Abstract 

Using scaled variables we are able to integrate an equation valid for isotropic 
and anisotropic Bianchi type I, V, IX models in Brans-Dicke (BD) theory. We 
analyze known and new solutions for these models in relation with the possi- 
(2JT)! bility that anisotropic models asymptotically isotropize, and/or possess infla- 

tionary properties. In particular, a new solution of curve (k ^ 0) Friedmann- 
Robertson- Walker (FRW) cosmologies in Brans-Dicke theory is analyzed. 
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1 INTRODUCTION 

The universe is nowadays at big scales homogeneous and isotropic as measured in 
the CMBR by the COBE satellite [|]], and must also has been having these proper- 
ties since, at least, the era of nucleosynthesis 0. In order to explain the isotropy 
of the universe from theoretical anisotropic models, many authors have considered 
the Bianchi models that can in principle evolve to a Friedmann-Robertson- Walker 
(FRW) cosmology. It has been shown that some Bianchi models in General Rela- 
tivity (GR) tend to their isotropic solutions, up some extent f|, and even that 
they can explain the level of anisotropy measured by COBE ||. Motived by these 
facts, we have been working in Brans-Dicke (BD) theory to investigate if Bianchi 
universes are able to isotropize as the universe evolves, and if its evolution can be 
inflationary. In previous investigations we have shown that anisotropic, Bianchi type 
I, V, and IX models tend to isotropize as models evolve f7|, ||. However, this may 
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happen for some restrictive values of u in the cases of Bianchi type I and IX, and 
only Bianchi type V model can accomplish an isotropization mechanism within BD 
current constraints M on u. It has been also shown that the isotropization mecha- 
nism in the Bianchi type V model can be inflationary, without the presence of any 
cosmological constant, when small values for the coupling constant uj are considered 



II , as in the case of some induced gravity (IG) models [[U], |TTJ. We have recently 
shown, however, that the isotropization mechanism can be attained with sufficient 
amount of e-folds of inflation, only for negative values of a;. In the present report 
we review and generalize some of the main results and present a new k ^ FRW 
solution in BD theory. 

This paper is organized as follows. In section |2| the BD Bianchi type I, V and 
IX equations are presented. In section |3] we review the main results on these models 
and present a new solution to curve k ^ FRW cosmologies. Finally, conclusions 
are in section |j. 

2 ANISOTROPIC EQUATIONS FOR BIANCHI 
MODELS 

In previous investigations 0, ^, [12|] we have used scaled variables, in terms of which 



our solutions have been given, therefore following we use them: the scaled field 
ip = (j>a 3 ( 1 ~ u \ a new cosmic time parameter drj = a~ 3u dt, ()' = the 'volume' 
a 3 = (1x0,20,3, and the Hubble parameters Hi = a/foi corresponding to the scale 
factors Oi = Oiirf) for i = 1,2,3. We assume comoving coordinates and a perfect 
fluid with barotropic equation of state, p = up, where v is a constant. Using these 
definitions one obtains the cosmological equations for Bianchi type I, V and IX 
models (in units with G = c = 1): 

(ijH t )' - ^R l3 = + (1 - u)u] po 3 ^ for z = 1, 2, 3. (1) 
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where a column sum is given by Rj = Sji?^-, where j =1, V or IX and 
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For the Bianchi type V model one has the additional constraint: H 2 + H 3 = 2H\, 
implying that a 2 and 0,3 are inverse proportional functions, note that 

the mean Hubble parameter, H = ^{Hi + H 2 + H 3 ), is for this Bianchi type model 
H = H 1 . 

Additionally, the continuity equation yields: pa 3( - 1+u ^ = const. = M v , M v being 
a dimensional constant depending on the fluid present. The vacuum case is attained 
when M v = 0. 

The system of ordinary differential equations, Eqs. (|TJ-^|), can be once integrated 
to getft 

^" " 3(1 ~ sh'-u) [m ^ (1 " 3u)v + 5 1 ^' + [2 + (1 - + 3z/V]m^ + 
_A_ [2 - 3zy + |(1 - vfuj] K(l - 3z/) r? + 5] 2 + (1 - v){h\ + h 2 2 + h 2 3 ) = , (5) 

where 5 is an integration constant, m v = §r|^, and the Hubble rates are written as 
follows (similar to the Bianchi type I model deduced in Ref. [|i~3|l ): 

h, = i ( g , + h 2 + g3 ) + ^ = ' - = (1 = 3 t"; " "w + 3(1 = ■ (6) 

3 3(1 — v)ip 

where the h^s are some unknown functions of r] that determine the anisotropic 
character of the solutions. If hi = for i = 1,2, 3 simultaneously, no anisotropy is 
present, which is the case of FRW cosmologies. Furthermore, Bianchi models obey 
the condition 

h! + h2 + h3 = Q (7) 

to demand consistency with Eq. (|6|). For the Bianchi type V model one has addi- 
tionally that hi = 0, since Hi = H as mentioned above. 

1 A similar equation, that is valid only for the Bianchi type V model, was derived in Ref. [ [l2| . 
Now we generalize that result in such a way that Eq. (^) is valid for Bianchi models I, V and IX, 
as well for FRW cosmologies. 



3 



Equations (Q) and (|7]) imply that the mean Hubble parameter is determined by 
ip alone: 



3H = H 1 + H 2 + H 3 



ip' (1 — 3u)m u r] + 5 
tp ip 



(8) 



In order to analyze the anisotropic character of the solutions, we consider the 
anisotropic shear, a = — (Hi — H 2 ) 2 — (H 2 — H 3 ) 2 — (H 3 — Hi) 2 , a = is a necessary 
condition to obtain a FRW cosmology since it implies Hi — H 2 — H$, cf. Ref. 
If the sum of the squared differences of the Hubble expansion rates tends 
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to zero, it would mean that the anisotropic scale factors tend to a single function of 
time which is, certainly, the scale factor of a FRW solution. 



The anisotropic shear becomes, using Eqs. (^|) and (|7|), 17(77) 
the dimensionless shear parameter using Eq. (Rf): 



3{hl+h%+hl) 
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If the above equations admit solutions such that a / H 2 — >• as 77 — > 00 (t —>■ 00), 
then one has time asymptotic isotropization solutions, similar to solutions found for 
the Bianchi models in GR [|J. In fact, one does not need to impose an asymptotic, 
infinity condition, but just that 77 ^> 77*, where 77* is yet some arbitrary value to 
warrant that a/H 2 can be bounded from above. 



3 ANISOTROPIC AND ISOTROPIC SOLUTIONS 



The problem to find solutions of Bianchi and FRW models in the BD theory has 
been reduced to solve the coupled system of equations @, (||) and ([!]). Let us 
present in the following subsections the known and new solutions. 

3.1 Bianchi type I model 

For this Bianchi model the known, the general solution is found in which the h^s 
are constants, then Eq. (H) is decoupled from Eqs. (H) and (J|), and the solution 
is ip = Aj if + Bj 77 + CV, where the constants Aj,Bj,Cj are reported elsewhere 
[0. [13| . This model can be solved in a general way since the curvature sum column 
Rj = 0, then Eq. (§) can be directly integrated. This solution represents a particular 
solution of Eq. (^). Direct substitution of ip into Eq. (|]) gives the Hubble rates, 
and into Eq. shows that solutions isotropize as time evolves, that is, a/H 2 — > 
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as 7] — > oo, see Ref. J7|. However, the isotropization mechanism is only possible for 
solutions such that A/ = B] — 4AjCi 

3(1 -i/) 2 (3 + 2o;) L 2(2 - 3i/) + 3(1 - vfuo v 7 

+ [2(2-3z/) + 3(l-z/) 2 c;]5 2 
+3(1 -u)\h 2 l + h 2 2 + hj)] (10) 

is negative 0. Then, some restrictions on u apply. For instance, in Dehnen's IG 
theory (Tl| u< 1, then the isotropization mechanism is not possible in this Bianchi 
model 101. 



3.2 Bianchi type V model 

This Bianchi model is more complicated because curvature terms are different than 
zero. Still, it is possible to find particular solutions of Eq. ([5]), since for this model 
the h^s are constants too, and this equation is decoupled as is the case of Bianchi 
type I model. The general solution of this model should be obtained through the 
general solution of Eq. (|]), yet unknown. We have found a particular solution that 
is again of the form ip = Ay rf + By r\ + Cy, where the constants Ay, By, Cy are 
reported in Refs. |7], [12| . This solution is such that 



-8(1 -3v 



:2 



18i/ + (l + 3z/) 2 a; 



K (11) 



is negative for uj > —18^/(1 + 3z/) 2 , so the solution tends to isotropic solution within 
BD theory constraints ||16|| , uj > 500, that is a/H 2 — > as rj — > oo. For this Bianchi 
type model, Dehnen's IG theory ]11| can achieve an isotropization mechanism ||. 



An inflationary behavior may be observed in type V models, but to get enough 
e-foldings of inflation (N ~ 68) one must demand that u < — § (12[ in consistency 
with previous results [ p]j . 



3.3 Bianchi type IX model 

Bianchi type IX model is the most complicated to solve, since curvature terms 
involve quartic polynomials of the scale factors, see Eq. @. For this Bianchi 
model it implies, by imposing the condition that h^s are constants, severe algebraic 
constraints on the scale factors, so it seems more likely that h^s are functions. This 
explains why no totally anisotropic (Hi ^ H2 ^ H%) solution has been found yet. 
In Ref. we have analyzed the case when the polynomial solution for ip is valid. 
In this case, unfortunately, we could not found explicitly the values of the constants 
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A IX , B ix ,Cjx- If this solution is valid, however, one has that h\ + h\ + h\ = 
Drj 2 + Ft] + G, where D, F, G are constants. Accordingly, Eq. (^|) indicates that 
the solution must tend, as time evolves, to the positive curvature FRW solution, i.e. 
one has again that a/H 2 — > as i] — > oo. However, a definitive answer will arrive 
by obtaining explicitly the values Ajx, Bi X , Cix- 

The only possible solutions for Bianchi type I and V models imply that the h^s 
are constants, whereas for type IX they are unknown functions of 77. An explanation 
of this fact resides in the property that Bianchi type I and V models have curvature 
terms of FRW type, whereas type IX has a very much complicated form, see Eq. @. 
So the things, it seems that the most general solution of Eq. fl5|) with h^s constants 
would give general solutions for Bianchi models I and V. The particular quadratic- 
polynomial solution of Eq. (||) represents in the case of Bianchi type I model its 
the general solution, whereas for Bianchi model V it is only a particular solution. 
Then, other particular solutions, possibly of non-polynomial nature, are expected to 
be found for Eq. fl£j|) that will reveal new aspects of Bianchi type V model. Finally, 
for Bianchi type IX model the quadratic-polynomial can be a possible solution, not 
yet confirmed. However, solutions with hi(rj) valid for the Bianchi type IX model 
are almost impossible to find because of the complexities involved in the curvature 
terms. 



3.4 FRW solutions 

It turns out that Eq. (|5]) is also valid for the FRW models when the anisotropic 
parameters vanish, i.e. hi — hi — hs — 0. Solutions of this equation solve FRW 
cosmologies in BD theory |7], 18]. The known solutions for ip are quadratic poly- 



nomials in i] as well. The general flat (k = 0) solution is a particular solution of 
Eq. (|5|) in which the coefficient of the quadratic polynomial term, A, is equal to 
the corresponding coefficient (Aj) of the Bianchi I case. For curved (k 7^ 0) FRW 
cosmologies the known particular solution is such that the coefficient of A of the 
quadratic polynomial term is equal to the corresponding coefficient (Ay) in the 
Bianchi V case. In this way, one can see a correspondence between anisotropic and 
isotropic solutions. 

We have found a new solution of Eq. @ that is valid for k 7^ FRW cosmologies. 
The new solution is: 



1 -3z/ 



2 



2 



^ = m " \T+3^) l^ + i 1 - 3 ") 1 !] |ki + 27/ + k 2 [«! + (1 — 3i/) 77] 1-3- J , (12) 

where ki,k 2 are arbitrary integration constants. This is the general solution of 
curved FRW cosmologies when the following relationships are valid: 

5 = 



6 



00 



-18 v 

:i + 3z/) 2 



(13) 



Though the latter relationship constrains the range of possible values of cu and u, 
one can find values of physical interest, e.g. u = — 1 that has some interest in string 
effective theories, see for instance Ref. ]nj. Moreover, when v — > —1/3 one obtains 
the GR limit u — > oo. Finally, one gets the dust model (f<just — 0) in the limit when 
uj — > 0, like in Dehnen's IG theory |TT| . 
The Hubble parameter is given by: 
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«1 + fEi + ^2 [Kl + (1 " 3l/) 17] 
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from which one can find the scale factor: 

Kx + 2-q + k 2 [ki + (1 - Zu) r]}^ 
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The BD field (0 = V>cr 3(1 ~ v) ) is obtained through Eqs. QT2p and (|TJ): 
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(16) 



Eqs. (15) and ([16]) imply that the sign of k 2 is equal to the sign of —k for 
most values of v. For open (k = —1) models this implies that k 2 must be positive, 
and for closed (k = +1 ) models k 2 must be negative which allows the solutions 
to (re)collapse: The value of « 2 determines the time of maximum expansion, so it 
is very related to the mass (m v ) of the model. On the other hand, K\ represents a 
?7-time shift. 

Because of the mathematical form of Eq. ( |14"D it is not possible to have an 
inflationary era that lasts for a sufficient time period to solve the horizon and flatness 
problems. To show the behavior of the models we have plotted the above-given 
formulae for different values of u, the curvature constant k, and integration constants 
k% and /%. We have chosen in all figures that |^=o= <P'o — as initial condition. 
Given a specific value for u implies two possible values of v consistent with Eq. (pf). 
For u — — 1, a value that makes BD theory to resemble string effective theories ||T9"f . 
it corresponds v\ — (2 — v^3)/3 ~ 0.0893 (that represents a quasi dust model) and 
u 2 = (2 + \/3)/3 « 1.2440. In figures fTJ and || we plotted the scale factor and BD 
field for V\. Figure [l] shows an open model with Ki = n 2 = —k = 1, whereas figure 
H is a closed model with K\ = k = 1 and k 2 = —0.001. In closed models, the smaller 
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k 2 is, the later in time they will recollapse. However, for models with 2/(1 — 3v) < 1 
one can find an upper limit on k 2 such that models never recollapse (this happens, 
for instance, when v = —0.3719, making u = 500, and if k 2 < —1); the effect of 
negative pressure avoids recollapse. 

One can compute the asymptotic limit for k = —1 models, when 77 — > 00. There 
are two limit cases, when the quantity 2/(1 — 3u) is greater or smaller than 1, cf. 
terms in the equations above. In the former case one has that 

^asimpl = {-k)^^ [(1 -3u)T]]^ 

fl — 3h>\ 2 3(1-./) 

0asim P i = m u K 2 ( j (-k) 2 ( 1 " 3 ") = const. , (17) 

whereas in the latter case one gets 

-2k 



6t-asimp2 
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2(l-3i/) -Zu l-9v 

(1 _ 3j/) (1-3./)* ^2(1-3^ 

1 — (1-3^ /— K 2 \ 2(l-3i/) 



2m,,- — - — ^——5 H77- r] 2 ^-'^) 2 . (1 



The plots of figure [1] tends to Eqs. (|T7[) ; one observes that tends to a constant 
value that can be fix to be G -1 through the right choice of the constants m„ and k 2 . 

Figures |3| and f| show an open and closed model, respectively, with uj = 500. 
Again, there are two possible values of u, u 3 = —0.2988 and z/ 4 = —0.3719. In these 
figures we have chosen the value of In figure ^ we have chosen K\ = k 2 = —k = 1, 
whereas in figure f| K\ — k — 1, k 2 — —0.8. The plots of figure |3] do not follow Eq. 
fll7l) nor fll8|) since 2/ (1 — 3v) is numerically very close to one, and the its asymptotic 
behavior lies somewhere in between. 

The relationship between the cosmic time t and rj, drj = a~ 3u dt, seems to be too 
complicated to be integrated in a closed form. However, in the asymptotic limits 
given by Eqs. fll?]) and fll8|) one obtains that 



i i 



tasimpl = (-k)^) (1~3U)— V — (19) 
-2 fc lra a _ 3zy) 2-( T ^ ; ) 2 2 _ gv(i+v) 



tasimp2 = 1 K2 J — T 77 2 ^) 2 . (20) 

p 2-9z/(l + z/) ' V ; 

In both cases the functions are monotonic, and for v\ and u 3 used in our plots, time 
grows as r) grows. Therefore, our time parametrization is appropriate. In this limit, 
one can express our solutions in 77 in terms of t to get that: 

^asimpl y/~—kt 

fl — 3u\ 2 3(1-/) 

0asim P i = ra v I — k 2 (-k) 2 ^~ 3 ^ = const. (21) 
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Figure 1: The scale factor and BD field a function of time 77. The plots represent 
an open model with k± — k 2 — —k = 1, and v = v 1 that implies u — — 1. 
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Figure 2: The scale factor and BD field a function of time 77. The plots represent a 
closed model with K\ — k — 1, K2 — —0.001, and v — v\ that implies uj — — 1. The 
model eventually recollapses. 
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Figure 3: The scale factor and BD field a function of time 77. The plots represent 
an open model with k± — k 2 — —k = 1, and v — u 3 that implies uj = 500. 
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Figure 4: The scale factor and BD field a function of time 77. The plots represent a 
closed model with k± — k — 1, K2 — —0.8, and v — v 3 that implies u = 500. The 
model eventually recollapses. 
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Eq. is a particular known flat solution in BD theory |2(J or the k = — 1 vacuum 
solution of GR. Eq. (^) is the flat space, Nariai solution 2~1] in which uo is given 
by Eq. 



4 CONCLUSIONS 

We have presented a set of differential equations written in rescaled variables that 
let us integrate a general equation to get Eq. (||), valid for Bianchi type I, V and 
IX models, as well as for FRW models. This equation is coupled to Eqs. ([!]) and 
(H), but for Bianchi type I and V models the anisotropic parameters (/ij's) are 
constants (since their curvature terms are of FRW type) and Eq. (|j) is decoupled. 
This property allows one to find the general solution of Bianchi type I model and 
a particular solution of type V; both solutions are quadratic polynomials. This 
solution let the models isotropize as time evolves, however, this can happens only 
for some parameter (u, v) range. The polynomial solution may also be valid for 
Bianchi type IX, but it is not proved yet. If it were, isotropization would be also 
guaranteed. 

We have found a new solution of Eq. (H) valid for curved (k ^ 0) FRW cos- 
mologies, that is, a solution with hi = 0. This is the general solution subject to 
the constraint given by Eq. (0). Accordingly, we have analyzed two cases of phys- 
ical interest: the case when u = —1, having some interest in string cosmology, 
that implies an equation of state of a quasi dust model (z/^0), and the case when 
uj = 500, consistent with current BD local experimental constraints [|l(J, implying 



that v ~ —1/3. Although in a different manner, both the scale factor and <fi grow 
as a function of the parametrized time 7] in all figures presented here. One can find 
specific values of the constants v, k 2 to have "closed" models without recollapse. For 
this to happen 2/(1 — 3^) must be less than 1, i.e., v must be negative. This is a 
known effect of negative pressures. 

The new solution is non-inflationary and for asymptotic times is of power-law 
type for the scale factor. On the one hand, when the quantity 2/(1 — 3z/) is greater 
than 1 the BD field tends to a constant value, then BD's dynamics evolves similar 
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to GR's. One the other hand, when 2/(1 — 3u) is smaller than 1, both scale factor 
and BD field behave asymptotically with a power-law, and the solution is equal to 
Nariai solution for flat space. 

Further solutions of Eq. (^) are in order, which can be either within FRW 
cosmologies or Bianchi type V or IX models. In particular, the general solution with 
Ttj's constants should provide the general solution of both curved FRW cosmologies 
and Bianchi type V model. 

Finally, our results could be also of interest for a;(0)-theories, where the value 
of the coupling parameter in some early cosmological era could have been rather 
different than its value nowadays, ui > 500. 

References 

[1] Bennett, C. L. et al. (1996). Astrophys. J. 464, LI; Gorski, K. M., et al. ibid. 
Lll; Hinshaw, G. ibid. L17. 

[2] Hawking, S. W. and Taylor, R. J. (1966). Nature (London) 209, 1278; Barrow, 
J. D. (1976). Mon. Not. R. Astron. Soc. 175 359. 

[3] Collins, C. B. and Hawking, S. W. (1973). Astrophys. J. 180, 317. 

[4] Barrow, J. D. and Sonoda, D. H. (1986). Phys. Rep. 139, 1. 

[5] Barrow, J. D. (1995). Phys. Rev. D 51, 3113. 

[6] Brans, C. and Dicke, R. (1961). Phys. Rev. 124, 925. 

[7] Chauvet, P. and Cervantes-Cota, J. L. (1995). Phys. Rev. D 52, 3416. 

[8] Mimoso, J. P. and Wands, D. (1995). Phys. Rev. D 52, 5612. 

[9] Cervantes-Cota, J. L. and Chauvet, P. A. (1999). Phys. Rev. D 59, 043501. 

[10] Fakir, R. and Unruh, W. G. (1990). Phys. Rev. D 41, 1783; ibid. 1792. 

[11] Dehnen, H. Frommert, H. and Ghaboussi, F. (1992). Int. J. Theo. Phys. 31 109; 
Dehnen, H. and Frommert, H. (1993). Int. J. Theo. Phys. 32 135; Cervantes- 
Cota, J. L. and Dehnen, H. (1995). Phys. Rev. D 51, 395; ibid. (1995). Nucl. 
Phys. B 442, 391. 

[12] Cervantes-Cota, J. L. (1999). Class. Quant. Grav. 16, 3903. 

[13] Ruban, V. A. and Finkelstein, A. M. (1975). Gen. Rel. Grav. 6, 601. 

[14] Chauvet, P., Cervantes-Cota, J. and Nunez- Yepez, H. N. (1991). In Proceedings 
of the 7th Latin American Symposium on General Relativity and Gravitation, 
SILARG VII, D'Olivo, J. C, et al., eds. (World Scientific, Singapore), p. 487. 

[15] Wainwright, J. and Ellis, G. F. R. (1997). Dynamical Systems in Cosmology 
(Cambridge University Press, Cambridge). 



14 



[16] Will, C. M. (1993). Theory and experiment in gravitational physics (Cambridge 
University Press, Cambridge) chapter 8. 

[17] Levin, J. J. (1995). Phys. Rev. D 51, 462; Levin, J. J. and Freese, K. (1993). 
ibid. D 47, 4282; Levin, J. J. and Freese, K. (1994). Nucl. Phys. B 421, 635. 

[18] Morganstern, R. E. (1971). Phys. Rev. D 4, 282; Ruban, V. A. and Finkelstein, 
A. M. (1976). Astrofizika 12, 371; Lorenz-Petzold, D. (1984). Astrophys. Space 
Sci. 98, 249; Barrow, J. D. (1993). Phys. Rev. D 47, 5329. 

[19] Dick, R. (1998). Gen. Rel. Grav. 30, 435. 

[20] Lorenz-Petzold, D. (1984). In Solutions of Einstein's equations: Techniques and 
results, (Springer Verlag, Berlin), p. 403. 

[21] Nariai, H. (1968). Prog. Theor. Phys. 40, 49. 



15 



